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ABSTRACT 

Poynting flux driven outflows from magnetized rotators are a plausible explanation 
for gamma-ray burst engines. We suggest a new possibility for how such outflows might 
transfer energy into radiating particles. We argue that the Poynting flux drives non- 
linearly unstable large amplitude electromagnetic waves (LAEMW) which "break" at 
radii r t ~ 10 14 cm where the MHD approximation becomes inapplicable. In the "foam- 
ing" (relativisticly reconnecting) regions formed during the wave breaks the random 
electric fields stochastically accelerate particles to ultrarelativistic energies which then 
radiate in turbulent electromagnetic fields. The typical energy of the emitted photons 
is a fraction of the fundamental Compton energy e ~ fhc/r e with / ~ 10~ 3 plus ad- 
ditional boosting due to the bulk motion of the medium. The emission properties are 
similar to synchrotron radiation, with a typical cooling time ~ 10 -4 sec. During the 
wave break, the plasma is also bulk accelerated in the outward radial direction and at 
larger radii can produce afterglows due to the interactions with external medium. The 
near equipartition fields required by afterglow models maybe due to magnetic field re- 
generation in the outflowing plasma (similarly to the field generation by LAEMW of 
laser-plasma interactions) and mixing with the upstream plasma. 



1 



1. Introduction 

A Poynting flux driven outflow from a magnetized 
rotator is a promising paradigm for gamma-ray burst 
(GRB) engines and there have been various imple- 
mentations of this concept (c.f. Usov 1992; Thomp- 
son 1994; Blackman et al. 1996; Meszaros & Rees 
1997; Kluzniak & Ruderman 1998). Such models are 
appealing in light of the fact that neutrino driven 
GRB emission from compact object mergers (Ruffcrt 
& Janka 1998; Janka et al. 1999) might fall short 
in radiative luminosity. Also, Poynting flux mod- 
els may provide a source of magnetic fields which 
could help alleviate the field amplification problem 
in GRBs themselves and in afterglows which requires 
fields of the order of equipartition (Wijers & Galama 
;i998)| , |Frail et al.(1999)j pranot et al.(1999)| ). 

Scenarios that could produce Poynting flux dom- 
inated outflows (PFDOs) require a source of mag- 
netic fields p> 10 15 Gauss, and rotation speeds of or- 
der ~ 10 4 /sec. The total available energy from a 
compact (Ro ~ 10 6 cm, M > Mq) object is then 
MVL 2 Rq/2 > 10 53 ergs and a dipole type luminos- 
ity 2B%R§n 4 /3c 3 - 2 x 10 50 ergs/sec. Such com- 
binations could be generated in a number of plau- 
sible cases: an accretion torus surrounding a black 
hole that formed from a neutron star merger (Rees 
& Meszaros 1997), strongly magnetized neutron stars 
("magnetars" Duncan & Thompson 1992) possibly 
formed from accretion induced collapse of a white 
dwarf with possible dynamo amplified fields in the hot 
young neutron star (Usov 1992; Duncan & Thomp- 
son 1992, Blackman et al. 1996), failed supernova 
lb (Woosley 1993), rapidly rotating magnetized black 
holes undergoing a Blandford-Znajek type energy ex- 
traction from supernovae (Lee et al. 1999, Brown et 
al. , 2000) or hyper-accreting accretion disks (Popham 
et al. 1999), with an MHD dynamo (Araya-Gochez 
2000). 

Without further specifying the nature of the cen- 
tral engine we suppose that GRB are due to such 
stellar mass objects, so that the conditions required 
for the generation of PFDOs are satisfied. Here we 
are interested in how Poynting flux models accelerate 
particles and produce gamma-ray emission. 

The suggestion that GRBs are PFDOs whose en- 
ergy is converted to particles and escaping radiation 
only at large distances from the central rotator po- 
tentially accounts for some important characteristics 
of GRBs (c.f. Blackman et al. 1996, Rees & Meszaros 



1997) (i) mass loading - PFDO can be launched 
without much matter, and particularly without much 
baryon contamination (as in pulsars, the wind may 
carry no baryons at all); (ii) compactness problem 

- the PFDO converts its energy effectively into par- 
ticles and gamma-rays at distances many orders of 
magnitude from the central engine, (iii) collimation 

- PFDO have a preferred direction along the rota- 
tion axis of the central engine; (iv) jets observed in 
pulsars and (probably) AGNs provide examples of op- 
eration of PFDO; (v) the energy from Poynting flux 
can be easily converted into high frequency electro- 
magnetic radiation. The presence of strong magnetic 
fields from the rotator provides electromagnetic en- 
ergy which might be tapped or converted into the 
fields which accelerate the radiating particles. 

PFDOs may differ from the conventional blast 
wave models in several important respects. The lepto- 
baryonic content of the outflow in the inner region is 
expected to be small by analogy to pulsar winds, and 
the energy carried by matter in the central region will 
be less than the energy carried by EM fields. As in 
the relativistic pulsar wind, the only other well estab- 
lished PFDO, the baryons may be absent altogether 
while energy flux in the electron-positron pairs may 
be as small as one millionth of the Poynting flux. 

The conventional estimates of the relativistic ex- 
pansion Lorentz factors and density of the outflowing 
plasma used in the framework of blast wave models 
may not be applicable to the Poynting dominated out- 
flows because of the strong magnetic fields (Paczyiiski 
1990, Usov 1999). Strong fields dramatically re- 
duce the importance of pair production from photon- 
photon interactions relative to the photon-magnetic 
field interactions. Other exotic QED processes, like 
photon splitting (e.g., Baring 1991, Adler 1971), will 
also become important for magnetic fields larger than 
the critical value of 4.4 13 Gauss. At this moment, we 
are unaware of the calculations of the Compton scat- 
tering cross section in supercritical magnetic fields, 
but we expect that it will also be strongly suppressed. 

The above suggests that in the central regions of 
PFDOs, the magnetic pair production and photon 
splitting may be more important than pair produc- 
tion in photon collisions. Note that this situation is 
realized in pulsars. There acceleration of particles to 
energies beyond 10 6 MeV proceeds relatively quietly, 
without producing a large number of observable pho- 
tons. Very few rotationally powered pulsars are ob- 
served at high energies. Thus conventional estimates 
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of the Lorentz factors and densities in GRBs based 
on two photon pair production may be irrelevant. 

To estimate a "revised" density of the generated 
pairs in the strong magnetic rotator case applica- 
ble to GRB, we use the analogy with pulsars where 
the number of pairs produced is related to the local 
Goldreich-Julian density ncj = VLB / {2irec) (Goldre- 
ich & Julian (1969)). The pair number density is 
usually A ~ 10 3 — 10 5 times larger than uqj in the 
inner regions. The factor A depends sensitively on the 
magnetic field, its curvature and accelerating poten- 
tial. A highly curved magnetic field produces denser 
and less relativistic outflows. 

In this paper we address the fundamental but un- 
resolved issue of how Poynting flux is converted into 
gamma-rays and the subsequent emission characteris- 
tics. In section 2 we give the basic picture of the mag- 
netized rotator magnetosphere and the propagation of 
large amplitude electromagnetic waves. In section 3 
we discuss the breaking of these waves, and the pre- 
dicted emission from random electromagnetic fields is 
discussed in section 4. There we address the spectrum 
and time dependence of the emission. The afterglow 
possibilities from our picture are addressed in section 
5. We conclude in section 6. 

2. Basic Picture 

Here we discuss the basic proposed scenario for the 
evolution of PFDOs. If there is a strong overall mag- 
netic dipole component, then near the central engine 
the PFDO will resemble a pulsar. Along the magnetic 
poles it will produce jet-like flow with helical mag- 
netic field, oscillating on scale length of (c/O) ~ 10 6 
cm. Initially Poynting energy flux dominates over the 
particle flux by a ratio ~ 10 6 - typical for pulsars. 

Inside the "light cylinder" r LC = (c/D.) = 3 x 10 6 
and the magnetic field falls off as r~ 3 , while beyond 
the light cylinder the magnetic field is dominated by 
the toroidal component which falls of as r -1 . Thus 
at large distances r 3> The we have 

B^B^{^) 2 (1) 

where Rq = 10 6 cm and Bq = 10 15 G are the size of 
the central engine and the initial magnetic field. 

Following the analogy to pulsar winds (e.g. Usov 
1992; 1994), at small distances the magnetic field of 
the wind is frozen in the plasma. As the plasma flows 
out, the plasma density decreases in proportion to r~ 2 



reaching a radius r t , where it becomes less than the 
critical charge density n t required for applicability of 
MHD (Goldreich & Julian 1969; Michel 1969, 1994; 
Coroniti 1990; Melatos & Melrose 1996). Locally the 
critical charge density n t is equal to the Goldreich- 
Julian density uqj. The latter decreases only as r -1 . 

In what follows, we normalize the real density of 
the PFDO n to the Goldreich-Julian density at the 
light cylinder riGj{rLc)i that is 

n = ArlG7 = A ^('My = 4xl0 15 cm- 3 (2) 
2nec \ c ) 

Since the plasma density satisfies n(r) = n{rLc){i'LC l r ) 2 , 
while naj{r) — nGj{rLc){ r LC l r ), the plasma density 
equals Goldreich-Julian density at 

rt = \r LC = 10 14 cm, (3) 

where we have chosen A = 3 x 10 7 (see also Usov 1994; 
Blackman & Yi 1998). The density and magnetic field 
at the transition radius are 

TH = n GJ (r t ) = 10 8 cm- 3 , B = 10 6 G. (4) 

The physical parameters near the transition region 
are listed in Table [|. 

At r > r t the the wind field is transformed into a 
Large Amplitude Electromagnetic Wave (LAEMW). 
The dimensionless strength parameter of the LAEMW, 
uq = eE / (tocJI), near the transition point is 

QR 3 B a g 

VV = o 10 ' ( 5 ) 



Though formally the plasma at transition point is overdense, 
u> p > f! a presence of a small "external" magnetic field with 
strength larger than 10 — 3 Gauss will allow wave propagation. 

Table 1: Physical parameters near the transition re- 
gion 



B, Gauss 


10 6 


ujb, rad/sec 


2 x 10 13 


;y o - — 


2 x 10 9 


tl = — , cm 


1.5 x 10" 


= i/2fiwfl, rad/sec 


6 x 10 8 


r s = — , cm 


50 




3 x 10~ 5 






LL. 


5 x 10 9 
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This large value implies that the electromagnetic 
waves can drive particles to superrelativistic energies 
in one period of oscillation. We suggest that gamma- 
rays will be emitted at the point where MHD breaks 
down due to the overturn instability of LAEMWs. 

The interaction of LAEMW with the plasma bears 
similarity to laser-plasma interactions. It is well 
known from laboratory experiments that interaction 
of strong electromagnetic radiation with plasma leads 
to a number of violent instabilities (e.g. Kruer 1988). 
The most important instabilities include resonant ex- 
citation of plasma waves when the frequency of the 
laser light matches the plasma frequency and non- 
resonant ponderomotive instabilities. The fact that 
GRB outflows may consist of purely electron-positron 
plasma might seem to play a stabilizing role, however 
pair plasmas are also susceptible to parametric and 
modulational instabilities (Lyutikov & Green 2000). 
Here we point out that in addition to parametric in- 
stabilities of EM waves, there is also an overturn in- 
stability associated with the non-linear evolution of 
the wave packet as it propagates through the plasma. 
We describe this instability in the next section. 

Once the EM wave overturns, it creates a broad 
spectrum of random EM fields. It is an electro- 
magnetic analog of the foam of the deep water sur- 
face waves. Overturn of the initial coherent wave 
will create random EM fields and currents which will 
be strongly dissipated. In general, there will be a 
nonzero component of electric field along magnetic 
field. We will call this state a relativisticly strong 
electro-magnetic turbulence. Since the end result is 
the dissipation of the magnetic energy, this regions 
may also be called a relativisticly strong reconnect- 
ing cites. Particles will be accelerated in this EM 
turbulence to relativistic energies (in the rest frame 
of the foam). As the accelerated relativistic particles 
move through the EM foam they generate the ob- 
served GRB emission. The emission is due to the ac- 
celeration in the random E-B fields. Emission would 
emanate from multiple breaking regions. The sepa- 
rate spikes in the GRB profiles would be due to sepa- 
rate "foaming regions" - separate regions of the wave 
break. 

In many respects the "foaming regions" resemble 
optical shocks observed in the laboratory when an 
optical pulse steepens as it propagates in an optically 
active medium (Agrawal 1989). A wave overturn cre- 
ates an electromagnetic shock. Interaction of such 
EM shocks will unavoidably lead to a creation of inter- 



nal plasma discontinuities with strong dissipative ef- 
fects. The "foaming regions" play the role of "internal 
shocks" (c.f. Piran 1999 for a review) of GRB theory. 
This term is usually used to mean "internal processes 
in the blastwave which generate bursts" not necessar- 
ily a specific mechanism, which we are providing here. 
Alternatively, the "foaming regions" maybe thought 
of as relativisticly reconnecting regions, where the en- 
ergy of magnetic field is converted into particles in a 
relativisticly turbulent manner. 

It is often assumed that some kind of Fermi scatter- 
ing process operates at the GRB relativistic shocks. 
In the shock Fermi process, particle scattering occurs 
across the shock between "scattering centers" of col- 
lections of e.g. Alfvcn waves. Technically speaking, 
shock Fermi acceleration is not the likely mechanism 
for particle acceleration for generic relativistic shocks 
because for upstream magnetic field angles greater 
than 1/r to the shock normal, where T is the bulk 
upstream Lorentz factor, the particles would have to 
travel faster than c to get back upstream along field 
lines. Thus unless the shock is strictly parallel (un- 
likely situation) then shock Fermi acceleration won't 
work. 

Hoshino et al. (1992) studied 1-D models of 
positron acceleration in relativistic perpendicular shocks 
(where the field is perpendicular to the shock normal) 
in ion-electron-positron plasmas. Their mechanism 
involves the absorption, by positrons, of gyro-phase 
bunched magnetosonic waves emitted by protons re- 
flected at the shock front. For <~ 10% ion fraction, 
non-thermal tails are seen in positrons. They dis- 
cuss limitations of their simulations with respect to 
their choice of electron-position number ratio, the fact 
that the simulations are 1-dimensional, and that they 
do not consider oblique shocks (shocks for which the 
fields are not exactly perpendicular to the shock nor- 
mal). They point out that electric fields may play 
more of a role in particle acceleration in the latter 
case. In general, they agree that more work is needed. 

We entertain the possibility that electric fields may 
be important for acceleration, but not necessarily only 
across a thin shock. Rather, we suggest that as an 
LAEMW overturns, the acceleration, even in pair 
plasmas, might be done directly by turbulent elec- 
tric fields (see also Lembege & Dawson 1987 for over- 
turn discussion). The thickness of the turbulent re- 
gion may be much larger than the thickness of a typ- 
ical shock. The overturn may force the region to be 
highly strongly turbulent over many gyro-radii. 
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The picture of the "foaming regions" bears some 
similarity to the "plasma turbulent reactor" (Tsy- 
tovich & Chikhalev , Norman & ter Haar 1975; Rees 
1967). In both cases there is a strong interaction of 
plasma kinetic turbulence with electromagnetic waves 
which may lead to a formation of stationary distribu- 
tion of particles. But the "plasma turbulent reac- 
tor" has been studied only in the limit of weak turbu- 
lence theory, when the energy density of the turbulent 
plasma fluctuations is much smaller than the total en- 
ergy density of plasma and when plasma is optically 
thick to transverse EM waves. The assumption of 
weak turbulence allowed for a closed self-consistent 
solution for the plasma turbulent spectrum and par- 
ticle distribution. On the other hand, the turbulence 
in the foaming regions of GRBs should be very strong, 
with an energy density much larger than that of the 
plasma. Detailed properties of such superstrong tur- 
bulence such as the stationary spectra, are not known. 
Under certain conditions we will show that the emis- 
sion from such turbulence does not depend in on its 
detailed properties (see section 4). We will not re- 
strict ourselves to a particular model of the turbu- 
lence, but we do discuss Phillips-type and relativistic 
hydrodynamics-type spectra in section 4.4 and 4.5. 

Each overturning region produces a burst of gamma- 
rays. The temporal behavior of such a burst might 
be most similar to a FRED (fast rise exponential 
decay) curve, since the fraction of LAEMW energy 
going into turbulence and particle acceleration likely 
occurs more rapidly than the subsequent radiation 
drain. The overall duration of a microburst is of the 
order of micro seconds (Eq. (p4])). To account for 
the energetics (e m b ~ 10 48 ergs per microburst with 
of order 10 4 microbursts as discussed later in section 
4) the overturning region should draw energy from 
a volume V ~ 10 36 cm 3 , i.e. the typical size of the 
turbulent region should be ~ 10 12 cm which is much 
smaller than the distance to the central region. Two 
other arguments may influence our estimates of the 
emitting region: (i) using the analogy with surface 
foam we can expect the transverse dimensions of the 
foaming region to be larger than the longitudinal di- 
mension, (ii) during a break of a wave packet energy 
from many wave periods is "piled up" in a narrow 
breaking region, reaching the thickness of one wave- 
length ~ 10 6 cm. This value also turns out to be of 
the order of c times the radiative decay length of a 
typical particle. 

The total number of overturning regions will de- 



pend on the structure of the flow, which is in turn de- 
termined by the structure of the central engine. For 
example, if the PFDO has a large dipole component, 
then the structure of the flow will be more regular and 
we expect only several regions of wave overturn. On 
the other hand if the structure of the magnetic field 
near the central engine is more multipolar, we expect 
many " foaming regions" . 

At the transition radius, momentum conservation 
dictates that there will remain some energy in the 
LAEMW which can be used for bulk plasma accel- 
eration. The bulk motion of the outflowing plasma 
should drive a relativistic shock wave into the am- 
bient gas. Synchrotron emission of particles at this 
shock could produce the observed afterglows. This 
would coincide with the conventional external shock 
models of afterglows (c.f. Piran 1999), while provid- 
ing a possible improvement: as we discuss further in 
Section 5, the outflowing plasma beyond the tran- 
sition radius may still be strongly magnetized, thus 
supplying the elusive magnetic field required in the 
synchrotron models of afterglows In this case the 
afterglows are something like the the analogs of the 
wisps observed around the Crab pulsar (the wisps are 
thought to be the location of the reverse shock where 
the pulsar winds pressure is balanced by the pressure 
of the supernova ejecta; Gallant & Arons 1994). Al- 
ternatively, the afterglows may be generated in a pro- 
cess similar to the one considered by Smolsky & Usov 
(2000), when the particles from the external medium 
are reflected from the strong magnetic field of the 
ejecta. In both cases however, we point out that the 
outflow can supply the magnetic field to the afterglow 
emitting region. 

3. Overturn of LAEMWs 

Here we consider the transition of the PFDO from 
MHD to the wave regime. In the MHD approxima- 
tion, the velocity of the wind equals the velocity of the 
field patterns. This wind velocity can be described 
as a wave phase velocity with trapped particles. In 
the laboratory frame this velocity is determined by 
the electric drift in the crossed electric and magnetic 
fields, v d = ^£ = f < 1. For the PFDOs, the 
leading EM terms are B^, Eq oc - so that the drift 
velocity (the velocity of the outflow) is directed along 

2 In the turbulent region magnetic fields with scales smaller than 
~ 10 6 cm are annihilated, leaving (or regenerating) larger scale 
fields which obey frozen-in condition. 
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the radial direction r. When the MHD approxima- 
tion breaks down, the phase velocity (i.e. the veloc- 
ity of the field pattern) will grow to v P h = Q > 1. 
Thus, during the transition from MHD approxima- 
tion to plasma wave solution the phase velocity of the 
field patterns will increase from subluminal to super- 
luminal. For plasma particles this process will look 
as if the EM wave is coming into plasma "from mi- 
nus infinity". As the EM wave penetrates into the 
plasma it exerts ponderomotive force along the direc- 
tion of its propagation which provides an acceleration 
for the plasma. Below we show that this process of 
acceleration is unstable toward wave overturn. 

The equations describing the evolution of the LAEMW 
waves in plasma are Maxwell's equations and the mass 
and momentum continuity equations for electron and 
positron fluids. To show the wave overturn simply, we 
consider a cold plasma. This is a reasonable approxi- 
mation when the particles acquire relativistic motion 
under the influence of the wave. 

Introducing the electromagnetic vector and scalar 
potentials A and <j), m the Coulomb gauge divA = 
the equations become (we use a system of units in 
which c = m = 1) 

(A - <9 t 2 ) A = 47r(n p v p - n e v e ) + d t V(f> 
Acf> = —4:ire(n p — n e ) 

dtPe,p + (V e ,p ' V)p e , p = 

±e (V</> + d t A - v 6iP x (V x A)) 



d t n e 



V (n e , p v e 



0. 



The force equations can be cast in the form 

dt (p e ,p T eA) - v e , p x (V x (p £:P ^ eA)) 
(±eV</> - V7) 



(6) 



(7) 



with the immediate consequence that the generalized 
vorticity, fl = V x (p e , P T eA) is conserved. If the 
plasma were initially quiescent without external field 
or vorticity, then condition (Q) guarantees that the 
vorticity vanishes at all times. Then @) simplifies to 



dt (p e ,p T eA) = (±eV0 - V7) . 



(8) 



Consider the propagation of a strong plane elec- 
tromagnetic wave along the z direction with a slowly 
varying amplitude a(£, z): 



a(i,z)e~ M ~ fe * z) 



(9) 



We assume that variations of the amplitude occur on 
the time scales much longer than the wave period 



and expand the dynamic equations in small quantities 
dta± and ?' a± . In nonlinear optics this is known as a 
slowly varying envelope approximation (e.g. Agrawal 
1989). To further simplify the treatment, we consider 
a circularly polarized EM wave for which the energy 
of particles in a wave remains constant during oscilla- 
tions (polar outflows from magnetized rotator should 
produce circularly polarized EM waves). For linearly 
polarized waves, an additional averaging over one pe- 
riod of oscillation is needed. In the zeroth order, we 
find that 



p_L ep =±eai. p z = const 



(10) 



Thus the transverse velocity of the plasma electrons 
and positrons are antiparallel and directed along the 
instantaneous magnetic field, but perpendicular to 
the electric field of the wave. The longitudinal ve- 
locity is arbitrary. 

To first order in the small parameters we find 



a z = 0, = 

dtPz + <9 2 7 = 0, 7 
d t n + d z (nv z ) = 



yJl + rf+Pl (11) 



(12) 



and a Schrodingcr type equation for the transverse 
vector potential (e.g. Berezhiani et al. 1992). Note 
that in the unmagnetized pair plasma, the variation 
of amplitude along the direction of wave propagation 



does not create charge separation, so n 



and 



Pz = Pz e>p - 

Equations ( |TT| ) and (|l2|) should be integrated along 
the group velocity characteristics dz/dt — v g . As- 
suming that all quantities depend on the variable 
£ = z — v g t, equations ( |TT| ) and (jl^) become 



~Vgd^Pz - 

df: In n = 



-da 

d^Vz 







(13) 
(14) 



For nonlinear waves, the group velocity v g depends 
on the Lorentz factor of the particle flow and on the 
local plasma density, v g (n, 7). The general solution 
of the system (|l^-pT|) for varying v g {n,^) seems to 
be untractable, so we first analytically illustrate the 
overturn for the case of constant group velocity and 
later resort to numerical integration for the fully non- 
linear case. 



Integrating ( |13| ) along the trajectories v g — const, 
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we find (Clemmow 1974) 



4. Emission in a random EM field 



n g (l + v 2 + y/1- (-l + 7i g 2 ) v A 

v z = — 

1 + Tig 2 + 1>Q 2 

7 _ n g 2 + ^l-(-l + n g 2 )u 2 (i5) 
n g 2 — 1 

where n g — l/v g and is the dimensionless am- 
plitude of the wave. Equation ([l5|) shows that for 
the intensity of the wave above the threshold value 

vo=v a l \Jl ~ v 2 (at which point v g = v z = — 7==) 

the Lorentz factor 7 becomes undetermined. This im- 
plies that hydrodynamic description of the medium is 
no longer valid: the waves overturn (Dawson 1959, 
Lembege & Dawson 1987; Lembege & Dawson 1989). 
During wave overturn, a mixing of various parts of the 
wave occurs, which destroys its oscillatory structure. 
Note that in the overturning regions the singularity 
in the density (E3) at this approximation is marks the 
regime of strong dissipation, where the bulk motion 
is randomized. It is thus a regime of particle acceler- 
ation. 

To further illustrate the wave overturn, we numer- 
ically integrate (O) and ( |l2] ) for a given profile of the 
incoming wave (Fig (Q)). 

The overturn of strong nonlinear waves is a well 
known phenomena in fiber optics where optical shocks 
(sometimes called kink solutions) form when a nar- 
row intense laser pulse propagating though a nonlin- 
ear medium (e.g. Agrawal 1989; Rothenberg & Gr- 
ishkowsky 1989: Trippenbach & Band 1998). Unlike 
relativistically strong LAEMW, the typical nonlinear 
intensities of the short laser pulses are small and usu- 
ally a nonlinear Schrodinger equation is used to con- 
sider wave breaking (Agrawal & Headly 1992). 

In an electron-ion plasma, the ponderomotive force 
felt by ions and electrons is different, resulting in dif- 
ferent accelerations and the creation of large unsta- 
ble currents. Generally electron-positron plasmas are 
more stable to parametric and modulational instabil- 
ities than ion electron plasmas since some of the non- 
linear 3rd order currents cancel out. However for a 
magnetized pair plasma the ponderomotive force act- 
ing on electrons and positrons is different. Thus if the 
outflowing plasma of GRBs has a large ion compo- 
nent or large nonoscillating magnetic fields, then the 
propagation of LAEMWs will be significantly more 
unstable. 



4.1. General remarks 

In this section we consider generation of X-ray and 
gamma-ray emission from the acceleration and propa- 
gation of relativistic electrons/positrons through rela- 
tivistically strong electromagnetic turbulence. Accel- 
eration of particles is due to random alignments of the 
electric field along the particle's velocity. Radiation, 
on the other hand, is due to the random components 
of the magnetic and electric field perpendicular to the 
particle's velocity. 

The typical electric field E ~ B produced near the 
transition radius is orders of magnitude larger than 
the Dreicer field in the mildly relativistic plasma with 
critical density n GJ : E D ~ ^ ~ l 1 "^) B. For elec- 
tric fields larger than Dreicer field, Coulomb collisions 
between electrons can be neglected and they are freely 
accelerated by the field. 

The spectrum of the plasma particles will be deter- 
mined by the competition between acceleration and 
emission in the stochastic EM fields. If these two pro- 
cess balance, a stationary spectrum will be reached. 
In this case the energy of particles would be radiation 
limited and only weakly dependent on the details of 
the acceleration mechanism. 

Acceleration of particles may proceed in two regimes 
which we refer to as "maximally efficient" and "stochas- 
tic acceleration" respectively. When the typical scale 
of the field fluctuations is larger than the typical ra- 
diative loss length, then the particle acceleration can 
be thought of as acceleration in an almost constant 
DC electric field and thus will be maximally efficient. 
If, on the other hand, the typical radiative length is 
much larger than the coherence length of the field fluc- 
tuation, then the particle will encounter multiple elec- 
tric field fluctuations and reversals before attaining 
steady state energy. Here the acceleration is stochas- 
tic. We expect that electromagnetic turbulence to be 
broadband, extending over scales from the skin depth 
up to the size of the turbulent region. The relative 
importance of the maximally efficient and stochastic 
acceleration types depends on the (unknown) spec- 
trum of the strongly relativistic turbulence. 

Emission from a relativistic particle in random E 
and B fields may be qualitatively related to emission 
from an electron in a random B field alone. To see 
this we note that the acceleration due to the electric 
field along the particle's velocity E || v is much less 
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(by a factor ~ j 2 ) than the acceleration due to trans- 
verse component of the electric field Ej_. (Landau & 
Lifshitz 1975). The component of the electric field Ej_ 
transverse to the particle's velocity acts similarly to a 
B field. Thus the emission of an electron in random 
E field will be very similar to the emission of an elec- 
tron in effective magnetic field of the same magnitude. 
Moreover, in a relativistic plasma the random electric 
and magnetic fields should be of the same order. 

The particle emission properties in turbulent mag- 
netic and electric fields also depend on the typical 
scales of field fluctuations. Relativistically moving 
particles emit most of the radiation along the parti- 
cle's velocity within an angle 80 ~ I/7. The resulting 
spectrum of the emission depends on whether the to- 
tal deviation angle of the particle during its transition 
through EM field is larger or smaller that ~ 1/7. For 
transverse acceleration, the relevant scale is the rela- 
tivistic Larmor radius tl = cy/u) e ff. Perturbations 
of the accelerating fields with scales larger than tl 
produce deflection angles larger than ~ 1/7. The re- 
sulting emission will be similar to synchrotron emis- 
sion. The total spectrum from an emitting plasma 
is then the average over the emitting volume of syn- 
chrotron spectra. 

On the other hand, for field fluctuations with wave- 
lengths smaller than the Larmor radius the deflection 
during emission is smaller than ~ 1/7. The particle 
moves almost straight along the line and experiences 
high-frequency jitter in the perpendicular direction 
from the random Lorentz force. Thus produces "jit- 
ter" radiation (Medvedev 2000, Landau & Lifshitz II) 
with a typical frequency ~ 7 2 c/A, where A is the scale 
length of perturbations. The spectrum is then deter- 
mined by random accelerations of the particle. As 
we show in Section 4.3 however, the expected turbu- 
lent scales are much larger than the Larmor radii, so 
that the emission will resemble synchrotron, not jit- 
ter, emission. 

4.2. Order of magnitude estimates 

To illustrate the some of the main points discussed 
above, we first consider the case when the smallest 
scale of the field fluctuations is assumed to be larger 
than the radiative loss length, which in turn is as- 
sumed to be larger than the non-relativistic Larmor 
radius. In this case, the acceleration and emission 
result from quasistatic electric and magnetic fields. 

Also in this case, the typical frequency and en- 



ergy loss of a particle follows from the analogy to 
synchrotron radiation in the an effective magnetic 
field B e j f = V< B 2 > + < E 2 >, where < B 2 > and 
< E 2 > are the rms of the fields: 



Ut~rUeff, (16) 

and where u e f / = e v / <- B2 >+< ga > j s the " effective cy- 
clotron frequency" . The total emissivity from a par- 
ticle is 



dU 
~dt 



2 2 



(17) 



In a steady state, the radiative losses would be ex- 
actly compensated by acceleration in the fluctuating 
electromagnetic fields. Assuming maximally effective 
acceleration, we find the energy gain to be 



dU/dt 



1 



V2 



mc v e ff, 



(18) 



where the factor 1/V2 is due to the fact that only 
the electric field contributes to acceleration. Equating 
Eqns. ([It]) and (|lj) we find the typical Lorentz factor 
of the particles to be 



(19) 



where 77 = c/uj e ff is the effective Larmor radius and 
r e — e 2 /mc 2 is the classical electron radius. The typ- 
ical emission energy is then given by a combination 
of fundamental constants only! 



-L—^ 40MeV. 
y/2 r e 



(20) 



Thus the typical energy of emission of a relativistic 
particle in turbulent EM field turns out to be inde- 
pendent of the spectral details of the turbulence in 
the maximally efficient regime. The energy of the 
particles is radiation limited. 

To see this more clearly, we compare the energy 
density in the particles in a stationary state with the 
equipartition energy density that would be achieved 
if approximately half of the energy of the initial mag- 
netic field were transformed into particle energy. The 
ratio of the magnetic to particle energy density may 
be expressed as a ratio of the relativistic Larmor and 
skin depths 



Urn 



1 



< 7 > \ r L 



» 1, (21) 
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( r ) 

where r s = y/< 7 >r, 
relativistic skin depth and the relativistic Larmor ra- 
dius and the latter relation follows for the relevant 
parameter choices. Since both quantities are proper- 
ties of the plasma as a whole, we used the averaged 
Lorentz factor of the particles < 7 > in their defi- 
nitions. Had we assumed equipartition, r\ ~ 1, the 
average energy < 7 > would have been 



(r) 

and r y L =< 7 > 77, are the 



< 7 >r 



10 9 



(22) 



which is much larger than ([H]). Thus the typical en- 
ergy density in particles is radiation limited and is 
less than the energy density in EM fields by a factor 
77 - 10 4 . 

The typical synchrotron decay time for a single 
particle is very short 



5 x 10" 8 sec 



(23) 



Yet in a steady state the energy lost by a particle 
due to radiation is rcsupplied by the energy stored 
in turbulent EM fields. The typical radiation decay 
times for the turbulence then will be 77 times longer 



f]t r 



5 x 10 4 sec 



(24) 



This radiation decay time maybe related to the dura- 
tion of micro spikes/bursts - each micro spike being 
due to one overturning region. 

4.3. Short turbulence scales: qualitative ac- 
count of stochastic acceleration. 

The estimates given above are the upper limits on 
the emission frequency and typical energy of particles. 
They are reached only if the correlation length of the 
field fluctuations l c is of the order of the radiative 
length 

77 r 2 

lr ~ c-f— = (25) 
dU/dt 7r e v ' 

which in the case of maximally efficient acceleration 
becomes l r rL^J^-- If, on the other hand, the correla- 
tion length of the fluctuating electric field is much less 
than the radiative length, l c <^l r: the effective accel- 
erating field will be lower by a factor yjl c /l r , which 
represents an effective RMS deficit in the coherence 
of the particle velocity and accelerating electric field. 



Suppose that the acceleration is lower than the 
maximum acceleration (Eq. (Ill) by a factor / < 1 



dU 

— ~ fu eff mc . 



(26) 



Then we would have a typical frequency of emission 



v t ~ fc/r e , 



a typical energy 



<7>= v/7-, 
and radiative length 



TL ( TL 



1/2 



(27) 



(28) 



(29) 



Suppose that the turbulent energy is concentrated 
near the smallest turbulent scale, which may be ap- 

(r) 

proximated by the relativistic skin depth = 
7 >c/lu p . (This will provide a lower limit to the 
efficiency of acceleration.) The ratio of the radiative 



length to the smallest scale of turbulence r 
becomes 



(r) 



» 



1 



73/2 rsre 



then 



(30) 



The efficiency of the acceleration / can be estimated 



Then and (BOJ) give 



r e TL 



2/7 



4 x 10" 6 . 



(31) 



Since at the transition radius ^5- = the efficiency 
of the acceleration / becomes 



/ = 



1/7 



2r P fl 



1/7 



2 x 10~ 3 . (32) 



The typical energy of emitted photons in this case 
c 



80keV 



(33) 



and the average energy < 7 >= \f]r^jr e = 5 x 10 3 . 

Relation (|33| ) gives the peak energy in the plasma 
frame. In the laboratory frame this should be mul- 
tiplied by the Lorentz factor of the flow. Since the 
observed peak energy are around 500 keV, we in- 
fer that the bulk Lorentz factors are relatively small 
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Tbuik ~ 10 (see also Stern 1999). | The low bulk 
Lorentz factors may reflect the fact that MHD accel- 
eration usually is not very effective (Michel 1969). 

We would like to stress that Eq. ([32]) was derived 
by assuming that the turbulent energy is concentrated 
on the skin depth scale. The presence of larger scale 
electric fields would increase the efficiency factor /, 
so ([52]) is a lower limit. The weak dependence |^| of 
this lower limit on the magnetized rotator's spin f2, 
its independence of the rotator's magnetic field and 
thus independence of the local plasma density param- 
eter, A, are interesting properties of our model. Said 
another way, when this lower limit on / applies, the 
energy of the emitted photons is independent of two 
out of three free parameters of the model (B and A), 
and only very weakly dependent on the third J7 1 / 7 . 

Next we confirm our earlier assumption that the 
emission resembles synchrotron emission not "jitter" 
emission. The typical nonrelativistic Larmor radius 
tl ~ 10 _3 cm is much less than the smallest turbulent 
scale (given by the the relativistic skin depth = 
c v / < 7 >/u p ~ 10 3 cm, < 7 >= 5 x 10 3 . Thus the 
particles will produce synchrotron-type radiation in a 
spatially varying magnetic field. 

We can also verify that that the plasma in optically 
thin to Thompson scattering (tt — OTnr t ~ 10~ 3 ) 
and to synchrotron selfabsorbtion at peak frequency 
(Pacholczyk Eq. 3.45). 

4.4. Spectrum of turbulence 

The actual spectrum of the relativisticly strong 
EM plasma turbulence has not yet be determined, so 
we resort to general arguments. One possibility is a 
power-law Kolmogorov-type spectrum in with dom- 
inant turbulent energy on large scales, transferred 
to small scales. The problem is that relativistically 
strong EM turbulence may not satisfy the following 
necessary conditions required for by Kolmogorov-type 
turbulence: (i) the turbulent cascades should be local 
in k space and (ii) there should be an inertial range in 
k space where resistive damping is weak. Both these 
conditions may not be satisfied: strong turbulence 
may be nonlocal and damping of turbulent energy by 
plasma particles maybe important. Q 

3 Recall that in our model the bulk Lorentz factors do not require 
Fbuik > 100. 

4 There are three free parameters:f2, Bo and multiplicity A 
5 Dettmann & Frankel (1996) argued that in relativistic hy- 
drodynamic turbulence far, where r is the scale parame- 



Another possibility is a Phillips-type spectrum 
(Phillips 1958). Following Phyllips (1958) we can ar- 
gue that at the point of overturn the profile of the 
wave becomes a discontinuous function (discontinuity 
of the zeroth order). The large wave number limit of 
the Fourier transform of such function will be oc fc _1 . 
Thus we can argue that the resulting spectrum of the 
EM turbulence will be £y oc fc~ 2 . 

Assuming that £u = Ck~ 2 and that the smallest 
scale of the turbulence is given by the relativistic skin 
depth, we find 

£ = J (PkE\\ - C J dk = Ck max = Cp (34) 

In this case, the maximum power per range of dk is 
oc k. Most of the turbulent energy is indeed con- 
centrated at the smallest turbulent scales ~ r s . For 
Kolmogorov-type turbulence, £11 oc /c~ n / 3 so that 
most of the turbulent energy is concentrated at the 
large scales. 

The turbulent plasma may also produce emission 
due to nonlinear conversion of the turbulent EM 
fields into escaping radiation. The typical frequen- 
cies will be of the order of the plasma frequency, 
ijj p ~ 10 9 sec _1 , falling into the radio wave band. An- 
other possibility is the production of radio waves by 
the overturning wave itself - an analogue of the super- 
continuum generation in optical shocks (e.g. Zozulya 
et al. 1999). 

4.5. Temporal evolution and spectra 

The temporal evolution of the emission from the 
foaming region depends on the onset and temporal 
behavior of the electromagnetic turbulence. In the 
initial stages of the wave breaking the energy of the 
initially coherent EM wave may be deposited at large 
scales, comparable to the size of the breaking region 
(~ 10 7 cm). After the initial energy injection, the 
energy cascades to smaller scales through turbulent 
interaction, reaching the smallest scale of the turbu- 
lence, the skin depth. 

This turbulent cascade may or may not reach a 
quasistationary state. The typical time for the devel- 

ter. if this is true, then Sk oc k~ 4 (compare with Kolmogorov 
£}, oc k~ 11 / 3 . If this scaling is true the energy is equally dis- 
tributed on all scales scales. 

' This is different from the breaking of surface waves on deep wa- 
ter, which form unstable cusps with a discontinuous derivative 
(discontinuity of the first order). In the case of surface waves 
the spectrum of the turbulence is ex fc -4 (Phillips spectrum). 
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opment of the turbulent cascade maybe as large as 
the overturn time of the largest turbulent eddies, ~ 1 
sec, if the turbulent cascade is local in k space, or 
as small as the typical microscopic interaction time, 
~ 1/ujp = sec, if the turbulent cascade is non- 

local in k space. The properties of the relativistically 
strong plasma turbulence have not been investigated 
yet, so we leave the question of the quasistationary 
state open and discuss observational features of both 
cases. 

If it takes longer than the plasma synchrotron cool- 
ing time for the turbulent energy to cascade to the 
smallest scales, then the quasistationary turbulent 
spectrum may not be reached. In this case we expect 
that the observed emission from such non-stationary 
turbulence will be almost time independent until (and 
if) the turbulence cascades to the radiation length 
scale. The reason is that initially the coherence size 
of the electric fields is much larger than the radiation 
length so the initial particle acceleration is maximally 
efficient. Particles are accelerated to roughly the same 
energy (since the energy is radiation, not acceleration, 
limited), so that the emission spectrum is quasista- 
tionary and similar to the synchrotron spectrum from 
a mono-energetic particle distribution. When the tur- 
bulent cascade reaches the radiation length, the accel- 
eration ceases to be maximally efficient. The particles 
are then accelerated in stochastic electric field, so that 
their spectrum will be power law with time dependent 
diffusion coefficient. The diffusion coefficient will be 
decreasing with time, acceleration will be less and less 
efficient and correspondingly the observed spectrum 
will be softening with time. 

In the limit of a very slow energy cascade, the radi- 
ation losses may drain the energy contained in the EM 
turbulence before the typical turbulent scale reaches 
the radiation length. In this case the spectrum of the 
particles and of the emitted radiation will resemble 
the spectra of synchrotron cooling sources. 

If on the other hand the turbulent energy cascades 
to small scales faster than the plasma synchrotron 
cooling time, a quasi stationary distribution of elec- 
trons will be reached. Acceleration due to small scale 
electric fields will be stochastic while acceleration due 
to large scale electric fields will be of DC type. The 
particle acceleration will depended on the spectrum 
of the EM turbulence. If most of the energy is concen- 
trated at large scales, as in Kolmogorov turbulence, 
then acceleration will be maximally efficient and the 
spectrum of particles will be approximately monoen- 



ergctic. If, on the other hand, most of the energy is 
concentrated at small scales, the acceleration will be 
stochastic. 

The observed spectra of GRBs will depend on the 
distribution function of emitting electrons, which, in 
turn, depends on several unknown parameters like 
the prevailing type of acceleration (stochastic or DC- 
type), proper boundary and initial conditions (contin- 
uos ejection of particles at low or high energies or "no 
flux" equilibrium states (Melrose 69, 71, 1980, Tade- 
maru 71, Park and Petrosyan, Katz 1994)), relative 
importance of escape of highly energetic particles in 
the long duration bursts. The stationary "no flux" 
equilibrium (Tademaru 71) will produce a relativis- 
tic Maxwellian distribution f(-/) oc j a e~' 1 ^°, with 
a > 2, while "reflecting boundaries" condition of Mel- 
rose (69, 71, 1980) will produce power law at 7 > 70. 
DC-type acceleration will tend to produce strongly 
peaked (monoenergetic) distribution, while stochas- 
tic acceleration will lead to spectral broadening. In 
addition, if particle escape is important, the result- 
ing spectra will be powerlaw-like, depending on the 
energy dependence of the escape probability (Melrose 
1980). [] This variety of possible spectra may serve 
well to explain the unusual variety of GRB's spectra, 
which are often powerlaws at small and large energies 
and sharply peaked near the spectral break (Band et 
al. 1993). We leave a more detailed consideration of 
spectral properties for a subsequent work. 

5. Regeneration of magnetic fields and after- 
glows 

The interface between the accelerated plasma and 
the ambient medium should lead to the formation of 
an external shock at large distances from the magne- 
tized rotator. This external shock could correspond 
to that inferred to be responsible for GRB afterglows 
(c.f. Meszaros & Rees 1993; Sari et al. 1998; Piran 



1999; Wijers fc Galama 1998;Frail et al.(1999), Gra 



not et al.(1999)), at large distances. To fit the af- 



terglows with the synchrotron emission of relativistic 
particles large magnetic fields, of the order of equipar- 
tition are required. Generation of such large magnetic 



fields have presented an unsolved problem (Gruzi- 
nov(1999")| ) in the blast wave model. 



7 In this respect we also note that the "plasma turbulent reac- 
tor" (refs**), which bears some similarity to our model prob- 
lem, produces a power-law spectrum p -3 , consistent with the 
observed spectral index —2. 
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In the PFDO model there is a possibility of re- 
generation of magnetic fields due to interaction of 
LAEMW with plasma. It has been known that inter- 
action of LAEMW in the form of lasers with plasma 
produces both large and small scale magnetic field 
(e.g., Stamper et al. 1975; 1978, Sudan 1993, Max 
1980). Most commonly, magnetic fields are generated 
on skin depth scales (Pegoraro et al. 1997). The typ- 
ical regenerated magnetic fields may be as large as 
initial magnetic field B ~ 10 6 G. It has been argued 
that such small scale quasistationary magnetic fields 
may explain the polarization observed in the after- 
glows (Medvedev 2000). Alternatively, the magnetic 
fields may be also subject to coalecence instability 
(Finn & Kaw 1977, Sakai et al. 199), in which ran- 
domly oriented magnetic islands merge forming large 
scale structures. 

There is also a possibility of generating large scale 
magnetic fields, with correlation length of the order of 
the size of inhomogeneities of the "of the laser beam" 
by inverse Faraday effect (Berezhiani et al. 1997), 
which may generate fields of the order 



B 



eR ± 



(35) 



where 70 ~ v§ is the typical Lorentz factor of particles 
accelerated by a EM pulse with a transverse dimen- 
sion R±, in our case the transverse size of the foaming 
region. 

We mention these possibilities here flag for fu- 
ture work, and to highlight the need to further inves- 
tigate possible analogies between field amplification 
in GRB and that of laser driven plasmas. The con- 
nection is that both plasma systems may be driven 
by LAEMW. 

6. Conclusion 

The two main points of this work may be summa- 
rized as follow: (i) Poynting flux dominated outflows, 
which so far has been invoked only in the context of 
central engines, possesses internal instabilities which 
may also explain radiation generation, (ii) emission 
from stochastically accelerated particles in turbulent 
electromagnetic fields, when electric fields are as im- 
portant as magnetic fields in the acceleration and ra- 
diation process, is a viable mechanism for the GRB 
emission. 

Our approach may also naturally explain the re- 
sult that the peak of GRB emission varies only over 



a small range of value from burst to burst. (Brain- 
erd 1994, 1997). For our maximally efficient mode 
of acceleration, the peak, when boosted into ob- 
server frame, becomes e max ~ /S,c/r e Fb u ik- If this 
maximally efficient regime were operating, the bulk 
gamma-factor Tbuik would not be much larger than 
10. The location of this peak is a very weak function 
of the parameters of the underlying rotator, / ~ Q 1 / 7 
and is independent of the progenitor's magnetic field. 

This can be compared to phcnomcnological inter- 
nal shock models in which the peak energy is propor- 
tional to -y^uJBTbuik, where j min is the low energy 
cutoff to the electron power law spectrum. All three 
quantities here arc usually taken to be free param- 
eters. It seems unlikely for their combination to be 
almost constant from burst to burst, and within each 
pulse of a given burst unless some physics dictates 
this to be the case. We have tried to add some of this 
physics. In our approach the remaining free parame- 
ter is the Lorentz factor of the bulk flow. 

PFDOs may also resolve the problem of the mag- 
netic field generation since in the Poynting flux driven 
outflows the large magnetic fields are supplied by the 
source. These fields do not necessarily have to be gen- 
erated in the external shock for the afterglow since the 
shock is a current sheet through which outflow and 
ambient particles mix (Smolsky & Usov 2000). At 
the same time, however there exists an unexplored 
analogy between field generation mechanisms in laser 
driven plasmas and GRB outflows that will have to 
be understood to determine the scale and structure 
of the outflow magnetic field. 

Other observational properties of the GRBs that 
may be explained in our framework: (i) GRBs show 
no correlation between the spectra and other mi- 
cropulse characteristics (e.g., Lee & Bloom 2000) - 
this is a direct consequence of the turbulent EM ac- 
celeration which produces photons with the frequency 
<~ c/r e ; (ii) the fact that pulses peak earlier at high 
frequencies and that bursts have shorter duration at 
higher energies may be due to the initial development 
of the turbulence: if the turbulence develops from 
large scales to small scales then initially the accelera- 
tion may be more effective since it is due to large scale 
electric fields (with coherence length larger than the 
radiative length). In this case initially particles are 
accelerated to the limiting energies 7 ~ \fT^Jr e while 
later, when the coherence length of the turbulence be- 
comes smaller than the radiative length, the particle 
distributions soften and radiation spectra emanates 
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at lower frequencies. The typical time for such a cas- 
cade should be of the order of the "vortex overturn 
time" on the largest scale of the turbulence, which 
may be on the order of seconds, (iii) the absence of 
correlation between the pulses and overall burst char- 
acteristics, interpreted as arising from random and 
independent emission episodes (Lee & Bloom 2000), 
is natural in our model since each wave overturn hap- 
pens independently, (iv) The temporal characteristics 
of the microbursts, FREDs, and the hard-to-soft spec- 
tral evolution is a consequence of synchrotron cooling 
of the reservoir of energy released during wave over- 
turn, (v) composite structure of GRBs (a burst being 
a sum of many independent emission events) (Stern 
& Svennson 1996) naturally follows from the model 
- each overturning region produces an independent 
microburst. (vi) the average power density spectrum 
of GRBs is well described as being due to selfsimilar 
turbulent-type process near marginal stability (Stern 
1999). This is reminiccnt of the cellular automata 
model of solar reconnecting regions (Lu & Hamilton 
1991) and may be related to the reconnecting "foam- 
ing" regions in our model. 

We also would like to point out that the particular 
mechanism of the wave instability, the wave overturn 
during MHD-wave transition, may not necessarily be 
the only one. Other plasma instabilities may con- 
tribute to the generation of EM turbulence. 

We would like to thank Michail Medvedev, Vladimir 
Usov, Norm Murray, Roger Blandford and Andrei 
Gruzinov for useful comments. 
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Fig. 1. — Profiles of (a) parallel momenta and 
(b) density in a wave. The profile of the wave is 
a = vq arctan(x — v g t) with v g — + fJ>n/j. Den- 

sity parameter /i = 0.1. Thin lines correspond to 
subcritical amplitude vq = 2.4, thick line correspond 
to ovcrcritial amplitude v = 2.5. 
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